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Abstract 

We present a supersymmctric model with a flavour symmetry 5*4 x Z3 and a CP 
symmetry which arc broken to Z^^^ in the charged Icpton sector and Z2 x CP in 
the neutrino one at leading order. This model implements an approach, capable of 
predicting lepton mixing angles and Dirac as well as Major ana phases in terms of 
one free parameter. This parameter, directly related to the size of the reactor mixing 
angle ^13, can be naturally of the correct order in our model. Atmospheric mixing 
is maximal, while sin^ ^12 ^1/3. All three phases are predicted: the Dirac phase is 
maximal, whereas the two Majorana phases are trivial. The neutrino mass matrix 
only contains three different real parameters at leading order and the neutrino masses 
effectively only depend on two of them. As a consequence, they have to be normally 
ordered and the model predicts the absolute neutrino mass scale and the sum of 
the neutrino masses. We study the main subleading corrections to the leading order 
results and show that they are in general small. 



1 Introduction 



The discovery of neutrino oscillations and the first information on lepton mixing angles 
has led to an intense research in this field. The data could shed light on the flavour 
mystery and help to establish some organizing principle relevant for fermion masses and 
mixing parameters. Before the determination of the reactor mixing angle ^13 several flavour 
symmetries, usually broken in a particular way, have been proposed, see [l], which give rise 
to mixing patterns with vanishing ^^13. A prime example is tri-bimaximal (TB) mixing [2] 
which was well compatible with the experimental data and could be derived from flavour 
groups such as ^4 ^ and 5*4 . The discovery of a relatively large reactor mixing angle [5] , 
however, requires most likely additional ingredients beyond the presence of a small discrete 
flavour symmetry. The simplest assumption which can be made is to invoke sufficiently 
large corrections to the leading order (LO) pattern such that current data [6] can be 
accommodated (?]. A disadvantage of this Ansatz is the fact that the predictability of 
such models is in general reduced. Discrete groups larger than A4 or 5*4 have also been 
considered in order to accommodate a non- vanishing ^13 [8j. A further approach consists 
in applying less stringent requirements on the residual symmetries in the charged lepton 
and neutrino sectors [o] . In this case 613 becomes essentially a free parameter. Correlations 
among this and the other mixing angles could still provide a testable prediction. 

The possible role of CP and its ability to allow for a non-vanishing reactor mixing 
angle, while constraining at the same time CP phases, has been explored to some extent 



in the past. In such cases CP acts non-trivially in flavour space 10 . A well-known 



example is given by the so-called fir reflection symmetry 11,12 , which exchanges a muon 
(tau) neutrino with a tau (muon) antineutrino in the charged lepton mass basis. If such a 
symmetry is imposed, the atmospheric mixing angle is predicted to be maximal, while ^13 
is in general non-vanishing for a maximal Dirac phase 6. Models combining 5*4 and CP 
can be found in the recent literature 13,14]. 



Taking these examples as starting point, we have recently considered an approach in 
which a discrete flavour and a CP symmetry are broken in a particular way such that 
the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix only depends on one free 
parameter. In the presence of a flavour symmetry Cf, CP transformations cannot be 
defined arbitrarily in flavour space, but are subject to certain consistency conditions [I8p9] . 
If these are fulfilled, we can consider a theory with Gf and CP and can study the case 
in which symmetries Ge and Gy are preserved in the charged lepton and neutrino sectors, 
respectively. Ge is a subgroup of Gf, while, generalising the idea of the /ir reflection 
symmetry, G^, is taken to be the direct product of CP and a Z2 symmetry contained in 
Gf. In 18 we have shown in a model-independent way that such a breaking leads to a 



mixing pattern with only one free parameter and thus mixing angles as well as CP phases 
are strongly correlated. We have furthermore shown that for Gf = 84, several new mixing 
patterns, well compatible with current experimental data |6], can be derived. 

In the present paper we discuss a concrete model in which the above idea of the break- 
ing of a flavour and a CP symmetry is implemented. The framework of the model is 
supersymmetric (SUSY) and the gauge group is the one of the Standard Model (SM). We 



^For other approaches comprising a discrete symmetry and CP see 15 - 17 . 
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focus on the lepton sector and consider as matter superfields only the leptonic ones of the 
Minimal Supersymmetric Standard Model (MSSM). There are additional supermultiplets, 
flavons and driving fields, which are singlets under the gauge group and are necessary 
for the correct spontaneous breaking of Gf and CP. We assume the mass of the scalar 
resonance found at LHC 20 ,21 can be accommodated by some mechanism. We consider 
G f = S4 X Z3 and a CP symmetry. The role of the symmetry is to facilitate the gener- 
ation of the charged lepton mass hierarchy without tuning. Two distinct sets of fiavons are 
present which contribute either to the charged lepton or to the neutrino mass matrix at LO. 
The vacuum expectation values (VEVs) of the former leave Ge invariant, while those of the 
latter G^. The vacuum alignment is achieved via vanishing of F-terms of the driving fields. 



The lepton mixing matrix is of the same form as in one of the cases exemplified in 18 
and the free parameter, directly related to the reactor mixing angle, is naturally of the 
correct order due to its origin from a subleading contribution to the neutrino mass matrix. 
The mixing angles agree within 3a with the latest global fit results [6] . Concerning the CP 
phases we find that the Dirac phase is maximal and the two Majorana phases vanish at LO. 
Since some of the contributions to the neutrino mass matrix which are compatible with 
the preservation of G^, are suppressed in our actual model, the neutrino mass spectrum 
can only be normally ordered. Indeed, the three neutrino masses effectively depend only 
on two different parameters so that, after fitting the two known mass squared differences, 
the absolute mass scale and the sum of the neutrino masses is strongly constrained. 



The paper is organized as follows. In section 2 we briefly recall the framework of 18 



and we apply it to the case of our interest. We assign the transformation properties of 
lepton and Higgs multiplets under Gj = x Z3 and CP and we determine the most 
general VEVs of the fiavon fields, compatible with the desired residual symmetries in the 
charged lepton and neutrino sectors. In section 3 we specify the superpotential of our 
model involving lepton fields, we analyze charged lepton and neutrino mass matrices, both 
at LO and at next-to-leading order (NLO), and we present the results for lepton masses 
and mixing parameters. In section 4 we discuss the vacuum alignment and we demonstrate 
that at LO the desired pattern of VEVs is a minimum of the scalar potential in the SUSY 
limit. Furthermore we examine the dominant NLO contributions to the fiavon VEVs. In 
section 5 we summarize the main results of our paper. Some details of the group theory of 
S'4 are contained in appendix A. 



2 Approach 



We elaborate on a proposal which we have recently presented in 18 . We assume a finite, 
discrete flavour group, here Gf = S4 x Z^, the direct product of the permutation group of 
four distinct objects and the cyclic symmetry Z3. This group is combined with a generalized 
CP transformation X which acts also on flavour space. If broken in a particular way, all 
lepton mixing angles and CP phases can be predicted in terms of only one free parameter. 

The group S'4 can be conveniently represented in terms of three generators S, T and 
U with S and U having order two, while T is of order three. For details about the group 
theory of S'4 see appendix |Xj A generic multiplet of fields is denoted by and transforms 
in a certain representation (r, q) of the flavour group Gf (r is among the S'4 representations 



2 



1, bJ or bj" 



with uj = e^^^^^). The CP transformation X 



and q stands for the charge: 
acts on as 

(j)'{x)=X<j)*{xcp) (1) 

with xcp = (x*^,— x). The choice of X is not arbitrary, but X has to fulfill certain 
conditions: first, we require X to be a unitary and symmetric matrix. 



XX^ = XX* = 1 



(2) 



so that CP^ = 1 holdsj^ Secondly, X also has to fulfill a consistency condition involving 
transformations which belong to the group Gf. For all elements g of the group Gf the 
following statement has to be valid: let A be the (unitary) matrix which represents g in 
(r,g), then an element g' E Gf should exist such that 



{x-^Axy = A' 



(3) 



with A' denoting the (unitary) matrix associated with g' in (r, q). Notice that g and g' are 
in general distinct elements. As has been shown in 18 , see also 19 , the mathematical 



structure of the group Gcp comprising Gf and CP is a semi-direct product and can be 
written as Gcp = Gf y\ Hcp with Hcp being the parity group generated by CP. 

Closely following the analysis which has been performed in 18 , we show that the lepton 
mixing parameters can be predicted from the breaking of the group Gcp to G^ = Z2 x CP 
in the neutrino and to Ge = in the charged lepton sector. Z^^ denotes here the 
diagonal subgroup of the group Z3 generated by T and the Z3 symmetry present m. Gf, 
but not contained in 5*4. |^ The requirement that Gy should be a direct product of Z2 and 
CP leads to a further constraint on X 



XZ* -ZX = 



(4) 



• q). 



with Z being the generator of the Z2 symmetry in the representation 

In 18 we have performed a comprehensive study of all admissible X for the possible 
choices of Z in case of Gf = S4. As can be checked these transformations X are also viable 
for G f = S4X Z3, since the abelian factor Z^ in Gf does not impose any further constraints 
on X. In order to simplify our consideration regarding the diagonal subgroup Z^^\ we 
choose a basis in which the generator T is diagonal in all representations of S^, its explicit 
realizations as well as those of the other two generators 5* and U can be found in appendix 

m 

By choosing Z = S, a consistent choice for X satisfying eqs.([2l [sl |4|) is 



Xi — 1 , Xi/ 



■1 , Xo 



X, 
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(5) 

Here we label with Xj. the realization of the CP transformation X in the representation 
r of S'4. Eqs. (T][5) specify the action of CP in flavour space. We observe that the 



18 



^For a possible generalization of this condition see 
^Compare also 23 for such a type of residual symmetry 
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transformations in eq. (|5| correspond to the choice X31 = 1 in [18], as it can be checked by 
changing basis from the one adopted here to the one chosen in |18| (via a unitary matrix). 
The generahzed CP transformation X3/ in eq. ([s]) is known as fir reflection symmetry in 



the hterature 11, 12 . It is interesting to notice that for an X satisfying eqs.([2| [3| |4j), also 
ZX does so. Therefore, also ZX provides an admissible CP transformation. 

We require that left-handed leptons / transform as (3', 1) under S4 x Z3. If the residual 
symmetry of the charged lepton sector Ge is Z^^^ , the combination m\mi [mi is the charged 
lepton mass matrix in the so-called right-left notation, I'^irti I) is diagonal and thus does not 
contribute to the lepton mixing, up to possible relabeling of the generations. For neutrinos 
we assume that they are Majorana particles and that their masses arise from the Weinberg 
operator. Preserving G,^ = Z2 x CP in the neutrino sector requires that the neutrino mass 
matrix rriy fulfills the invariance conditions 

Z^m^Z = , Xs/m^Xg/ = m* . (6) 



It follows that 
with 



UpMNS "^v UpMNS = mf"'^ , (7) 




2cos^ V2 2sme 

UpMNS = I — cos6' + i-\/3 sin^^ -\/2 — sin 6* — 2-^3 cos | K 
- cos sm — sm cos 6 

and mf"'^ containing the neutrino masses rrii, i = 1,2,3. The unitary matrix K is diagonal 
with entries ±1 and ±i which encode the CP parity of the neutrino states. Notice that, at 
this level, the mixing matrix is only determined up to permutations of rows and columns, 
since neither charged lepton nor neutrino masses are constrained by the above conditions. 
All mixing parameters, angles and Dirac as well as Majorana phases, are predicted in terms 
of the parameter 6, which depends on the entries of the neutrino mass matrix and can take 
in general values between and tt. From eq.(|8]) the mixing angles can be read off 

2 11 

sin^ 9i3 = - sin^ 9 , sin^ 6*12 = 7— — , sin^ 6*23 = - ■ (9) 

3 2 + cos 20 2 

Thus, for 6 ^ 0.185 the best fit value for the reactor mixing angle can be accommodated 
well and also the solar and atmospheric mixing angles are within their 3 a ranges [6] . For 



a X analysis we refer to 18 . Since the second column of Upmns in eq.feh is tri- maximal 



the solar mixing angle is constrained by sin ^^12 > 1/3 22 . The Majorana phases a and 
P are trivial, 

sina = , sin/3 = , (10) 
while the Dirac phase 6 is maximal 

|sin(5| = l (11) 
and thus \ Jcp\ = I sm29\/{6\/3). A characteristic feature which has already been noticed 



in the literature 11 12 of the mixing matrix in eq.rt8| is the fact that the absolute values 



^For the convention we use for mixing angles and CP phases see appendix in 18 



4 



of the entries of the second and third rows are equal, i.e. |?7^j| = \Uri\ for i = 1,2, 3. This 
directly leads to maximal atmospheric mixing, since sin 623 = cos 623 and for non- vanishing 
613 also to a maximal Dirac phase, cos 5 = 0. 

In the subsequent section we want to construct an explicit model and thus have to 
specify several aspects. First of all, we consider an extension of the MSSM in which the 
MSSM Higgs doublets h^^d are not involved in the breaking of the flavour symmetry, but 
gauge singlets only charged under Gf are responsible for this breaking. Secondly, we choose 
the right-handed charged leptons to transform as e*^ ~ (1', 1), fi'^ ~ and t'^ ~ 

under 5*4 x Z3. This choice does not affect the result for lepton mixing, but turns out to be 
useful for accommodating the hierarchy among charged lepton masses correctly. Indeed, 
as we argue below, see section [3} this hierarchy can be naturally achieved with a judicious 
choice of flavons. The transformation properties of the matter superfields and the MSSM 
Higgs doublets under 5*4 x Z3 and the auxiliary symmetry Ziq can be found in table [1} 
Notice that the role of this auxiliary symmetry is mainly to sufficiently disentangle the two 
flavour symmetry breaking sectors. 





I 
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hd 


5*4 


3' 


1' 
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1 


1 


1 


Z3 


1 


1 


u 




1 


1 


Zl6 


1 


ul 


1 


1 




1 



Table 1: Matter superfields and MSSM Higgs doublets h^^d of our model. As one can see, 
only one of the latter is charged under the auxiliary symmetry. Note that the auxiliary 
symmetry Ziq is effectively a Zg symmetry at the level of operators with matter superfields, 
since {IhuY carries a phase The phases are u = e^'^*/^, cos = e^'^*/^ and coie = e^'^^l^^ , 



Since we wish to break the group Gcp spontaneously to = Z./^^ and G^ = Z2X CP, 
it is convenient to list the most general form of the VEVs that fiavon fields in the different 
representations of 6*4 can take and which leave Ge and Gj, invariant, respectively. In order 
to preserve the group Ge = Z^^^ flavons ^ and ^' in one- dimensional representations 1 and 
1' of 5*4 having a non- vanishing VEV should transform trivially under Z3. In the case of 
a fiavon x being a doublet under 5*4 the group Ge can be preserved only if x carries a 
non-trivial charge under Z3. The VEVs are 

X ~ (2,u;) : (x) oc J ^ and x ~ {2,u') : (x) cx J ^ . (12) 

Flavons ip ^ 3 and ~ 3' can always acquire non-zero VEV, independent of their charge 
under Z3. The latter charge, however, determines which of the components of the fields 
actually acquires a non-vanishing VEV: 



^-(3,1) : (V^)(x ,V^~(3,u;) : (^) (x 1 ,z^~(3,u;2) - (^) ^ q (13) 
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and the same being true for the structure of the VEVs of (/? ~ 3' under S4. 
The VEVs that preserve G^, are 

(0 = , (O = ' (X) = ( ' (^) = I 1 I , (^) = I 1 I , (14) 

with the parameters and v^f, being real, v^i and being imaginary, while is in general 
complex. Since = Z2 x CP, it is irrelevant for the vacuum structure preserving Gi, 
how these fields transform under Z3. At LO only the flavons ^, x ^"^^ f can couple to 
the the Weinberg operator, see appendix |A} Their couplings are suppressed by A which is 
the generic high-energy cutoff in our model. Since it is also the scale which is relevant for 
the Weinberg operator it is in general inversely proportional to the light neutrino masses. 
Upon breaking to Gi, the most general form that the neutrino mass matrix can take is 



with y^, y^ and y^ indicating the Yukawa couplings to the flavons ^, x and ip, respectively. 
These are real because of the invariance of the original theory under a. CP symmetry 
generated by X. As expected, this matrix is diagonalized by Upmns in eq.([8]). The 
imaginary part of is responsible for the breaking of the fir symmetry of the matrix m^, 
and thus its relative size compared to the other VEVs is directly related to the size of the 
reactor mixing angle ^13. For this reason, this imaginary part arises effectively from an 
operator with two flavons in the model presented in the subsequent section, whereas the 
flavons transforming as 1 and 3' contribute via operators with one fiavon to the matrix miy. 
Notice that the vanishing of the imaginary part of promotes Gi, to Z2 x Z2 x CP. As 



one can check, the neutrino mass matrix in eq.(15) contains effectively four independent 
real parameters. This tells us that we cannot make a prediction for the neutrino mass 
spectrum in the general case. In our specific model, however, the contribution leading to 
an equivalent effect as the real part of is suppressed and the neutrino mass spectrum 
can only to be normally ordered. 



3 Model with ^4 and CP 

In this section we define our model. As anticipated, we adopt a SUSY framework, neglecting 
all SUSY breaking effects. The gauge group is that of the SM and the lepton fields are 
those of the MSSM. The flavour symmetry is G/ = 6*4 x and the CP transformation 
X is defined as in eq.([5]) for each representation (r, g) of the flavour group. We use as 
auxiliary symmetry a Ziq group in order to forbid unwanted operators. The transformation 
properties under GfX Z\^ of the leptons and MSSM Higgs doublets h^^ are collected in 
table 1, while those of the flavons and the driving fields are presented in tables 2 and 3, 



respectively. We discuss the results for lepton masses and mixing at LO in subsection 3.1 



and the those at NLO in subsection 3.2 All operators directly contributing to the lepton 
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mass matrices are non-renormalizable in our model and are suppressed by appropriate 
powers of the cutoff scale A. We consider two different sets of flavons one responsible for 
the breaking to Ge and one for preserving Gp. At LO the fields whose VEVs leave Ge 
invariant only contribute to charged lepton masses, while those leaving Gy intact generate 
neutrino masses. At the subleading level this changes and all flavons contribute to charged 
lepton as well as neutrino mass matrices. 



3.1 Leading order results 

In order to implement the approach described in section [2] we let the three generations 
of left-handed leptons U transform in the irreducible three-dimensional representation 3', 
while right-handed charged leptons e'^, mentioned above, are singlets under the 

flavour group. In this way, the generation of the hierarchy among the charged lepton masses 
is facilitated. At the same time, the three generations of charged leptons are distinguished 
by their transformation properties under the cyclic group Z3. In our model, charged lepton 
masses arise at the one-, two- and higher-fiavon level and thus are naturally hierarchical. 
We include the fiavon fields (pE ~ (3',w) and xe ~ (2,^) under S4 x Z3. The most 
relevant operators for charged lepton masses are ther 



We = yr{l'PE)r^hd/K + 1/^,1 (/<^|)/i^/irf/A2 + y^,,2{lXEVE)f^''hd/ A"^ ■ (16) 

Since we impose CP as a symmetry on the theory in the unbroken phase, all couplings in 
We and in the following superpotentials have to be real. As explained in section |4| at LO 
the flavons xe and ipE acquire the VEVs 

(Xe) = J ^ and (pe) = t;^^ [ 1 j , (17) 

which leave the subgroup Ge = ^3^'* of S4 x Z3 invariant. We find mi to be diagonal in 
flavour space with only non-zero entries in the (22) and (33) elements which correspond to 
the muon and tau lepton masses 




m 



V V 

(2 Vi^^iv^^ + yf,,2VxE) -jf- (hd) and = (ha) 



respectively. Their relative hierarchy is correctly reproduced for (note the VEVs of xe and 
ifE can be in general complex) 

I^xbU^^bI ~ with A ^0.2. (19) 

Due to the size of |f<^g| we also predict tan/3 = (hu) / (hd) to be of order 10. The electron 
remains massless at this level and its mass is generated via higher dimensional operators, 
discussed in section 13.21 



■^Adding the field xe is mainly necessary in order to achieve the correct vacuum alignment, as seen in 
section |4] 

^We mdicate the contraction to a trivial singlet 1 of 5*4 by (• • • ). 
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Xe 




C,N,C,N 


Xn 
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3' 
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3' 
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^3 
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1 
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1 
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■2^16 


1 


1 






Us 





Table 2: Flavons and their transformation properties under Gf x Ziq. The flavons labelled 
with E are charged under Z3 and neutral under Zig, whereas those with an index N only 
carry non-trivial charge under Zig. The phases are oj = e^'^*/^ and Us = e^'^'^l^ . 



As explained, neutrino masses arise in our model from the Weinberg operator. The 
fiavon multiplets ^at, and (f>^ are in (1, 1) and in (3', 1) under 5*4 x Z3, respectively, and 
couple at the level of one fiavon insertion, while the fields ~ (1', 1) and xn ~ (2, 1) 
couple together at the level of two flavons to neutrinos. Such a situation can be arranged 
by an appropriate choice of the charges under the auxiliary symmetry Zig, see table |2} The 
lowest dimensional operators responsible for neutrino masses are thus 

= y,,i{ll)hl^N/A' + y.AllK^N/A'' + y.,2{llvNK/A^ + y.AllXN^NK/A' ■ (20) 
As shown in section kl at LO the flavons S,N) ^n-, ^'ni Xn and v^at develop the VEVs 



{In) = , {^n) = v^^ , (^^) = v^:^ , (xn) = ^Xiv I J . {^^n) =v^^ ^1 j (21) 

with v^j^, v^^ and being complex, while v^'^ is imaginary. We also find that v^^^, v^^ 
and v^^ have a common phase (up to relative signs). The neutrino mass matrix rrii, takes 
the form 

1X1/ 

mi,= I -Uy-ixy 2uu + ixu ty - I ^""l' (22) 





with parameters 



tv = yu,l ^ ' = ^^^.2 -J^ , Xy = l ?/^,3 ■ (23) 



The common phase of f , v-^j^ and v^j^ is unphysical as regards lepton masses and mixing 
at this level, since it can be factored out from m^. Thus, the parameters t^, Uy and x^ 
can be regarded as real. Among the leading contributions coming from the operators in 



eq.(20) it is the one comprising ^'^ and Xn which renders the neutrino mass matrix in 



eq.(22) invariant only under Gy = Z2 x CP and not Z2 x Z2 x CP. The matrix nxy in eq. 



(22) is slightly different from the one in eq. (15) since, after extracting the overall phase, 
the combination Xn^'n only develops an imaginary VEV. The latter is responsible for a 
non- vanishing value of 9i^. This mixing angle is naturally smaller than the other two ones, 
because it arises from an operator with two flavons instead of only one as the other two 



8 



angles. In order to achieve the correct size of ^13 the relative strength of the contribution 



from the last operator in eq.(20) must be suppressed by a factor A with respect to the 
contributions from the first three operators in eq.(20). This is the case, if we choos^ 



AA . 



(24) 



As mentioned in section |2| the neutrino mass matrix rriiy is diagonalized by a mixing 
matrix of the form given in eq.(|8|. At the LO such a matrix is the PMNS mixing matrix, 
since there is no contribution from the charged lepton sector. The angle 6 is directly related 
to 6'i3, see eq.([9]), and depends on the two parameters and 



tan 26' 



Xy 



The neutrino masses read 

mi = \ti, + Uu\ 
where we have defined 



1712 



A ' A 

sign(Mi, cos2^^) a/3u^ + x1 



A 



(25) 



(26) 



(27) 



In order to show that neutrino masses as in eq.(26) imply normal mass hierarchy, we 
compute the atmospheric 



Ami 



— m 



{huf 
A 



(2^ 



and the solar mass squared differences 



mo 



A 



Since the solar mass squared difference is positive, we need 

Uy {2ty+Uy) < , 



(29) 



(30) 



while Am^^-j^ can be either positive (normally ordering) or negative (inverted ordering). An 



inverted hierarchy requires tyUy > which is incompatible with the requirement in eq. ( 30 ) . 



In contrast, a normally ordered spectrum entails tyiiy < 0. The two possible solutions 



{ty < 0, Uy > OT ty > 0, Uy < O) arc indeed compatible with the constraint in eq.(30), if 
we require 2\ty\ > \uy\. As one can see from eq.(29), we have to tune 2ty + Uy to be small 



in order to achieve a small value of Aml^J Aml^^.^. Since the neutrino masses only depend 
on two real parameters ty and Uy, both of them are fixed by the measured values of the 
mass squared differences and as a consequence the absolute scale of neutrino masses and 

^Choosing all flavon VEVs associated with the neutrino sector to be of order A is only convenient, but 
a priori not necessary. 
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also their sum are fixed. Taking the best fit values for Arrig^j and Am^^^^^ as found in [H] 

we get 

mi ^ 0.016 eV, m2 ^ 0.018 eV, mg^ 0.052 eV (31) 

and 

Em^^ 0.086 eV. (32) 

If we assume {hu) = 100 GeV and A = 10^^ GeV such that {hu)'^/A = O.leV, the param- 
eters ty, Ui, and Xi, are of order A and A^, respectively, for all couplings being of order 
one. 



3.2 Next-to-leading order results 



In this subsection we address the most relevant contributions arising from shifts in the 
LO vacuum and from operators with more flavons than those discussed above. For the 
following discussion of the subleading contributions it is important to keep in mind that 
we have chosen different expansion parameters in the charged lepton and neutrino sectors, 



see eqs. (19) and (24). We first discuss the most important corrections to the neutrino 



mass matrix. These arise from the shift in the VEV of the flavon xn, see eq.(54), if plugged 



into the forth operator of eq.(20), and lead to a contribution to the neutrino mass matrix 
in eq.(22) which is suppressed by A^ with respect to those coming from the first three 



terms in eq.(20). This new contribution cannot be absorbed via a re-definition of the three 
parameters t^, and x^, in rriy, but instead requires a forth parameter, called Pi, in the 
following. This parameter is real, since the correction to the VEV of xn as well as are 
both imaginary and all couplings are real. Furthermore, the contribution is /xr-symmetric 
in flavour space. Thus, the corrected neutrino mass matrix reads 



m 



NLO 



with 



ty + 2Uy -Uy " IX y + Py -Uy + iX y + Py 

-Uy - iXy + Py 2Uy + IXy +Py ty - Uy 

-Uy + iXy +Py ty - Uy 2Uy " IX y + Py 



Pu = I yu,3 a A 



{huT 
A 



(33) 



(34) 



being real, since i v^i^ is real. This matrix has the most general form which is compatible 
with the preservation of the group Gy in the neutrino sector (modulo the overall phase 
which also shows up in the parameter py and thus can be factored out). Consequently, 
the results for the mixing parameters remain unchanged, whereas the neutrino masses, see 



eq.(26), acquire corrections proportional to py. Since py is suppressed by A^ with respect 
to ty and Uy, the correction to the LO results for the neutrino masses is also of the order 
A2. 

Further corrections at order A^ arise, if we consider a shift in the vacuum of the field 
and re-compute the contribution coming from the forth operator in eq.(20). Corrections 



at order A^ appear if we include operators with fiavons out of which three belong to the 
set {xe,'^e} and one is either ^n, or (p^ and compute their contribution using the LO 
vacuum alignment of these fiavons. In both cases the structure of the neutrino mass matrix 
j^NLO jg ah^gj-gfj in such a way that it is no longer compatible with the preservation of the 
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residual symmetry G^. Thus, such corrections will also affect our results for the mixing 
parameters. However, these corrections are suppressed by at least with respect to the 



those coming from the operators in eq.(20) and can be safely neglected for this reason. 

In the same fashion we can analyze the subleading contributions to the charged lepton 
mass matrix which is diagonal at LO and does not allow for a non-zero electron mass. The 
latter instead is generated in two ways: first through operators with five flavons with an 
index N , if we take into account the shift in the VEV of the flavon xn (in the following we 
only mention those that actually give a non-zero contribution to the electron mass, omit 
all (real) Yukawa couplings and do not specify whether these operators arise from only one 
independent 5*4 contraction or whether there are several independent ones) 

le'hd^ljXN^N/A^ + le^hd^NXN^l^/A-^ + le^hd^NXk^N / (35) 

and secondly through operators with six flavons of this type (one of them being an 
even number of fields xn and the rest either ^at, or (^n), if we use the LO vacuum of 
these fields, 

+ le''hdeN^'NxlvN|^•'' + le'^h^j^xlvl/^' + le''hdi'^x%VN|^•'' ■ (36) 

All these operators induce contributions to the charged lepton mass matrix which are of 
the order {hd) in the {li) elements, i = 1,2,3. Notice that these are all equal. In this 
way, an electron mass of order X^{hd) can be achieved which is in agreement with the 
observed mass hierarchy among the charged leptons. The reason for these operators to 
give rise to a contribution which is equal for all three elements of the first row of mi can 
be easily understand by realizing that the form of the contribution can be simply derived 
from symmetry considerations, namely by studying the most general form of the VEV of 



a field in the representation 3 of which leaves the group G^, invariant, see eq.(14). 

Also the other elements of the charged lepton mass matrix are corrected at a subleading 
level. One type of corrections originates from the shifts in the VEVs of the fields xe and 



ipE, see eq.(56), and results in (31) and (32) elements of order A^ and (21) and (23) elements 
of order A^ in units of (hd) ■ A second source of corrections are operators with more flavons 
which turn out to have the following structure: an operator with n'^ comprises two fields 
from the set {xe, V^b} in order to saturate the Z3 charge and either one field xn together 
with three flavons of the type {C,n,^n,'^n} or three fields xn and one field ^n, C,n or (p^, 
while an operator with r'^ only requires one field from the set {xej'Pe}, while otherwise 
maintaing the same form as the operator with /i*^, e.g. 

llfUdxlil^NXN/h!" + Ir'^hdipEilxN/M' . (37) 

Taking into account the LO form of the VEVs of fields with index and we see that 
such operators contribute at the same level as the leading operators which are, however, 
evaluated with the shifted vacuum of xe and ipE- Eventually we find the order of mag- 
nitudes of the different entries of the charged lepton mass matrix including subleading 
corrections 

A^ A^ A6 

mfLo ^ I A^ AM {hd) . (3^ 
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The charged lepton masses follow the correct hierarchy 

me : : rrir X'^ : \^ : I (39) 

with tjIt- ~ (hd) and we can estimate the size of the contribution to the lepton mixing 
angles coming from the charged lepton sector as 

~ with i] = 12, 13, 23 . (40) 

Thus, such corrections have a small impact on the results for the lepton mixing parameters 
and can be neglected in our analysis. 

4 Flavon superpotential 

In order to construct the fiavon superpotential responsible for the alignment of the fiavon 
VEVs we assume the existence of an i?-symmetry U{1)r [25] under which matter superfields 
carry charge +1, fields acquiring a non-vanishing vacuum (i.e. hu^d and all fiavons) are 
uncharged and so-called driving fields, indicated with the superscript "0" carry charge +2. 
In this way all terms in the superpotential either contain two matter superfields (like the 
terms leading to lepton masses) or one driving field, hke those to be discussed in this 
section. In the limit of unbroken supersymmetry, the F-terms of the latter have to vanish 
and in this way the vacuum of the fiavons gets aligned. The driving fields necessary for 
our construction of the fiavon superpotential are given in table [3] Their transformation 
properties under the flavour and auxiliary symmetries are chosen in such a way that at the 
renormalizable level the vacuum of the sets of fields {xe, '^e}, {^n, ^n, Xn, Vn} and are 
aligned separately. In this way, we can enforce the first set of fields to spontaneously break 
Gcp to Zf^\ the second one to Z2 x and to G^. Since also the fiavon superpotential 
is invariant under CP, all couplings are real. 





SB ) SB 




SAT 


X% 




SAT 


5*4 


1 


2 


1 


2 


3' 


1 




UJ 


(jj 


1 


1 


1 


1 


Zie 


1 


1 


W4 


ojI 


oj'l 


1 



Table 3: Driving fields and their transformation properties under G/ x Ziq. The fields 
labelled by E are charged under Z3 and neutral under Zig. The opposite holds for fields 
labelled by A^. The phases are u = e^'^*/^ and co^ = e^'^*/^. 



Thus we can divide the fiavon superpotential Wfi in three parts Wfi^e, wji^^, and Wfi^^ 

Wfl = Wfi,e + Wfi^^ + Wf^ (41) 

with 

Wfi,e = de ^UxeXe) + ae ^U^e^e) + K (XeXeXe) + Cg (Xe^e^e) , (42) 
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+du (XnXnXn) + ey {x%^n^n) + fu iN{^%^N) + 9u {v%VnVn) (43) 

and 

Wfi,^ = l%M^ + ami'^i'^). (44) 

A few things are noteworthy: the driving fields and are defined in such a way that 
the former only couples to {xeXe) and the latter only to {lpe'^e)] similarly, we use the 
possibility to define the fields and in such a way that only couples to 937V at the 
renormahzable level; furthermore, notice that we neglect all terms containing the MSSM 
Higgs doublets hu,d, since these are completely irrelevant in the discussion of the alignment 
of the fiavon VEVs which typically are of order 10^^ GcV. 

For the F-terms of the driving fields with an index E we find 



dw 



fi 



dw 

dw 



-- 2aeXE,iXE,2 (45a) 

-- de(ip%^i + 2ipE,2'PE,3) (45b) 



be Xe,1 + Ce {<Pe,Z + '^VE,1<PE,2) (45c) 
be Xe,2 + Ce (<^B,2 + '^'P E ,l(p E ,3) (45d) 



dx%^ 



dx%,2 

Equating these to zero and solving for the VEVs of the fiavons xe and ipE get 



{Xe) ^v^e 5 ) and ((^e) = ^ 1 j 



(46) 



up to symmetry transformations belonging to S/^. Obviously, also the trivial vacuum is a 
solution which, however, can be easily excluded by requiring that at least one of the fiavons 
Xe and kpE acquires a non- vanishing VEV. The two VEVs and v<^^ are related through 



2 2 I An\ 



■^e 



with v^j^ being a free parameter which can also be complex. The presence of this parameter 
indicates a fiat direction which exists in the potential in the SUSY limit. Since v^^^ and v^j^ 
arc directly related it is natural to expect them to be of the order of magnitude, as we have 
done in the preceding discussion. Being their size related to a fiat direction we however 
cannot fix their absolute size, but choose this by hand in order to correctly accommodate 
the mass of the tau lepton. 
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Analogously, the F-terms of the driving fields Xn ip% are 



Cat + du^N + ^N^N + 2 &^ Xn,iXn,2 + c,, + 2</?iv,2</?Ar,3) (48a) 

0^1^ Xat,! + Gz. (¥^^,3 + '^Vn,iVn,2) (48b) 

c^i/ XAr,2 + {^%,2 + 2(^jv,i(^jv,3) (48c) 

fu iN^N,l + 2 ((^^_i - (PN,2^N,3) (48d) 

= /i. iN^N,3 + 2 9u {^lf,2 - fN,l<PN,3) (48e) 

= fu In<Pn,2 + '2gu {Vn,3 - ^n,iVn,2) (48f) 



9Wfl _ ^ ,,2 , „ t,^2 

9Wfl _ ^ ,^2 , ^ ^„2 
%^,2 

9Wfl _ f ^ ,^ , ^,„2 

dwfi 

9Vn,2 
dwfi 



Equating also these F-terms to zero we find several possible solutions for the vacuum. 
These can all be classified and have different characteristics: one class is formed by the 
trivial vacuum, one class by the vanishing of the VEVs of the fields xn and v^at, a third 
class comprises the vacua in which the VEV of xn vanishes, while the other fiavons acquire 
non-zero VEVs, and the last class requires {^n) — and all other VEVs to be different 
from zero. We are indeed interested in this last class and as one sees it can be selected by 
requiring that also the fiavon xn being in the two-dimensional representation of S\ acquires 
a non-zero VEV. Then all solutions belonging to this class are related by 6*4 transformations 
or by a relative sign between the VEVs of the two components of the fiavon xn- The latter 
statement can be easily understood, since we notice that Wji only depends quadratically 
on Xn- Given these different possibilities we choose as vacuum 



(49) 



and {^n) = 0. The VEVs v^^^, v^^ and v^^j^ are related through 

with v^^ parametrizing a flat direction in our fiavon superpotential. Notice that this 
parameter is in general complex so that also the VEVs of C,n, Xn and (/^^v are in general 
complex. Apparently, then CP should be broken in the neutrino sector. However, notice 
the following: for appropriately chosen parameters of the potential, the VEVs have all the 
same phase, up to tt. Thus, in the neutrino mass matrix this phase turns out to be an 
overall phase which can be factored out, as discussed above. Notice furthermore that (at 
least) at the LO level we can re-define the fields of the superpotential in such way that 
the latter acquires an overall phase which has no physical meaning. The fact that f^^, 
v^j^ and v^j^ are directly related through the parameters of the superpotential naturally 
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explains that they have the same order of magnitude which we choose to be A A. Eventually 
requesting the F-term of to vanish leads to 



2 



4 = . (51) 

For M'^/a^ > we immediately see that the VEV of is imaginary, as it should be in 
order to preserve G^. Notice for M ~ AA also {C,n) is of this order of magnitude. 

We now turn to the discussion of the non-renormalizable contributions to the fiavon su- 
perpotential. These induce shifts in the vacuum obtained at LO. Taking into consideration 
the sizes of the LO VEVs (fields with an index E have a VEV of order A^ A, whereas those 
with an index N have VEVs of order A A), we find the most relevant subleading corrections 
to arise from the following three operators 

Su,i {x%^'nXn)/A + Su,i {x%^'nXn)/A + Su,2 {(P%^'nXn^n) / ^ ■ (52) 

Adding these to w ji and recomputing the F-terms in a linear expansion in the shifts in the 
vacuum, we find that they induce a correction to the VEV of xn 

with Sv^j^- denoting the shifts in the i^^ component of xn, {XN,i) = + Sv^^^.. Being 
proportional to v^^ it carries the phase associated with the fiat direction and shared with 
the LO VEVs of the fields ^n, Xn and ipN and at the same time it is proportional to v^'^ 
which is imaginary. Its relative suppression with respect to the LO VEVs of the flavons 
with index is A. Thus, the shifted VEV of xn can be parameterized as 



with a real, showing that this shifted VEV, as expected, leaves invariant the group G 



see eq.(14). The other VEVs do not acquire non-trivial shifts at this order. Also this is 
expected on symmetry grounds. As one can check, all terms, which lead to contributions of 
order A^ A^ or larger to the F-terms of the driving fields, if LO fiavon VEVs are plugged in, 
only contain fiavons with an index N and thus do not introduce further non-trivial shifts 
in the vacuum. The most relevant terms including the flavons xe and which indeed do 
change the alignment are 

eKv'D/A + lli^lxE)/^ + d(xl)/A . (55) 
Their intrinsic size is A^ A^ and thus they are relatively suppressed by A*^ compared to the 



LO terms in eq.(44). Since the VEVs of xe and carry an undetermined phase such 
contributions lead to a shift in the phase of the VEV of so that it acquire a small real 
part which is suppressed by A^ with respect to f^j^. We neglect such a small correction in 
the discussion of the lepton mass matrices. 

Similarly, we can discuss the subleading contributions to Wfi^e- The ones leading shifts in 
the vacuum of xe and which cannot be absorbed into its LO form comprise six flavons. 
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Two of them belong to the {xe, V^b} and are necessary in order to saturate the symmetry, 
while the other four ones carry an index N (there are either three fields of the {C,n, ^n, ^n} 
and one power of xn or three fields xn and one of the set {^n,C,n,'^n})- Assuming the 
sizes of the VEVs, as discussed above, these contributions to the F-terms of the driving 
fields and x'e ^f the order A^A^ and thus suppressed by A"^ with respect to the 



contribution from the renormalizable terms in Wfi^e, see eq.(42). Consequently, they induce 
shifts in the VEVs of xe and at the relative order A^ such that the VEVs take the form 

(Xe) = j and {^e) = v^^ ^1 + -fX^ j (56) 

with /3, 7, e and rj having absolute values of order one and being in general complex. 



5 Summary 

In this paper we have constructed a SUSY model for leptons with a flavour group Gf = 
5*4 X Z3 and a generalized CP symmetry. The spontaneous breaking of these symmetries to 
a Z3 group in the charged lepton and to Z2 x CP in the neutrino sector leads to a mixing 
matrix with one free parameter. The latter is determined by the elements of the neutrino 
mass matrix and is in our model naturally small. This parameter is directly related to the 
reactor mixing angle which is predicted to be of the correct order A. At the same time, 
the other mixing parameters are also predicted, especially the atmospheric mixing angle is 
maximal as well as the Dirac phase. The two Majorana phases are trivial. The sine square 
of the solar mixing angle is slightly larger than 1/3. We predict the three neutrino masses 
to effectively depend on two parameters only at leading order and thus the mass spectrum 
can only be normally ordered with the absolute neutrino mass scale being fixed. We 
have discussed in detail the vacuum alignment of the flavour and CP symmetry breaking 
fields and have shown that this alignment is stable against the inclusion of higher-order 
terms in the flavon superpotential. We have computed the leading as well as subleading 
contributions to the charged lepton and neutrino mass matrices and have shown that all 
our predictions are only slightly corrected by the latter contributions. The most relevant 
effect of the subleading corrections is the generation of the electron mass. The charged 
lepton mass hierarchy is correctly reproduced without any tuning, since the different lepton 
masses arise from operators with one, two and several flavons, respectively. 



Note Added: During the completion of this paper another paper 26 has been sub- 
mitted to the arXiv that contains a similar SUSY model in which the flavour group S4 and 
a generalized CP symmetry are broken to subgroups Ge and G^ = Z2X CP in the charged 
lepton and neutrino sector. 
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A Group theory of S. 
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We follow the convention and notation of 24 and we repeat for the reader's convenience 



the form of the 5*4 generators S, T and U in this basis, the Kronecker products and the 
form of the Glebsch Gordan coefficients. 

The generators S, T and U in the five different representations are chosen as follows 

1:5 = 1, T = l, U = 

l':5 = l, T = l, U = 



1 


^ 









1 . 




-1 


2 


( 


2 


-1 




2 


2 




-1 


2 


( 


2 


-1 




2 


2 



S=l{ 2 -1 2 \ ,T=\Ouj^O\,U 





3' : S=l\ 2 -1 2 l,T=|Ow2ol,[/ 

with u! = e^'^*/^. These fulfill the relations 

52 = 1 , = 1 , f/2 = 1 , 

{STf = 1 , {SUf = 1 , {TUf = 1 , {STUy = 1 . 

In order to perform explicit computations of invariants under S4, we need the Kronecker 
products as well as Glebsch Gordan coefficients. The Kronecker products are of the form 
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1 X n = ^ y H , l'xl' = l, l'x2 = 2 
1' X 3 = 3' , 1' X 3' = 3 , 
2x2 = 1 + 1' + 2, 2x3 = 2x3' = 3 + 3', 
3x3 = 3' X 3' = 1 + 2 + 3 + 3', 3x3' = l' + 2 + 3 + 3'. 

In the following we list the Glebsch- Gordan coefficients using the notation a; ~ 1, x' ~ 1', 

(z/1,1/2)*, (^1,^2)* ~ 2, {zi,z2,z3y,{zi,z2,s3y 3, (2;, 4, 4)*, (5^,4,4)* ~ 3'. 

For a singlet multiplied with a doublet or a triplet we find 

iW X 2 : {xyi, xys)* ~ 2 , {x'yi, -x'ys)* ~ 2 , 

1^'^ X 3 : {xzi,xz2,xz3y 3 , {x'zi,x'z2,x'z3y 3' , 

iW X 3' : {xz[,xz'2, xz'^y ~ 3' , {x'z[, x'4, x'z'^y ~ 3 . 
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For a doublet coupled to a doublet 

2x2: yiy2 + 1/2^1 ~ 1 , Vih - 1/2^1 ~ 1' , (1/2^2, ViVif ~ 2 
and for a doublet multiplied with a triplet we have 

2X3: (z/l2;2+?/2^3, 1/1^3+1/2^1, 1/1^1+1/2^2)* ~ 3 , (l/l2:2-?/2^3,l/l^3-?/2^1,Z/l^l-?/222)* ~ 3' , 

and 

2x3': (?/i4-y24' Z/i4-?/2^i> yi^'i-1/24)* ~ 3 , (?/i4+l/24> l/i4+l/2^i> yi4+y24)* ~ 3' 
For the products 3x3 and 3' x 3' (with Zi to be replaced by 4) we get 

ZiZi + 2:253 + ^3^2 ~ 1 , (^1^3 + ^2^2 + Z^Zi,ZiZ2 + ^1 + 2:353)* ~ 2 , 

(2:2^3 — 2:3^2, 2:1^2 — 2:2^1, 2:3^1 — 2:1^3)* ~ 3 , 

{2zizi - 2:2^3 - 2:3^2, 2z3Z3 - 2:1^2 - 2:2^1, 2^252 - ^1^3 - 2:3^1)* ~ 3' , 
while the Clebsch Gordan coefficients for the product 3x3' read 

Ziz[ + 2:22:3 + ^3^2 ~ 1' , (2:12:3 + Z2Z2 + Z:iz[, -{Z1Z2 + 2:2^1 + ^3^))* ~ 2 , 

(22:1 z'^ — 2:22:3 — ^32:3, 22:3^3 — 2:12:2 — 2:22:^, 2^24 — 2:1^3 — 2:32:^)* ~ 3 , 
(2:22:3 — Z^Z^i Ziz'2 — Z2z\^ z^z'^ — Ziz'^Y ~ 3' . 
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